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Abstract 

Ph. 

We compute the damping rate of a fermion in a dense relativistic plasma at 
zero temperature. Just above the Fermi sea, the damping rate is dominated by 

x 

5_i ■ the exchange of soft magnetic photons (or gluons in QCD) and is proportional 

to (E — fj,), where E is the fermion energy and fi the chemical potential. We 
also compute the contribution of soft electric photons and of hard photons. 
As in the nonrelativistic case, the contribution of longitudinal photons is 
proportional to (E — n) 2 , and is thus non leading in the relativistic case. 
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The properties of quasiparticles in ultrarelativistic (UR) plasmas have attracted much 
attention in a recent past Q. A crucial property of a quasiparticle is its decay (or damping) 
rate: a quasiparticle which propagates in a plasma is not stable, as it undergoes collisions 
with the other particles of the plasma, and the very concept of a quasiparticle makes sense 
only if its damping rate is small enough. 

The damping rate of an electron propagating in a nonrelativistic (NR) plasma was com- 
puted almost forty years ago by Quinn and Ferrell ||, At first sight, this damping rate is 
infinite, due to the singular behavior of the Rutherford cross-section at small angles. How- 
ever Quinn and Ferrell were able to obtain a finite result because the Coulomb interaction 
is screened in a plasma (Debye screening), and in the case of a degenerate plasma, they 
showed that the damping rate is proportional to (e p — 6f) 2 , where e p = p 2 /2m is the NR 
kinetic energy and ep the Fermi energy, when e p is slightly larger than Bp. The damping 
rate remains finite for a non zero temperature T: only the value of the Debye screening 
length is modified. 

It is interesting to extend the calculation of the damping rate to the case of relativistic 
plasmas: one may have in mind either electromagnetic (QED) plasmas, such as found in 
white dwarves or in the core of nascent neutron stars, or chromodynamic (QCD) plasmas 
such as the quark-gluon plasma which is believed to be formed for large enough values of the 
temperature T and/or the chemical potential fi. The NR results are not easily transposed 
to the relativistic case, because the exchange of magnetic (or transverse) photons in QED 
or of magnetic gluons in QCD becomes important, while in the NR case it is suppressed 
by powers of (f/c) 2 with respect to the exchange of electric (or longitudinal) gauge bosons, 
and is usually neglected. These magnetic photons, or gluons, give rise to severe infrared 
(IR) divergences which are not easily cured, because there is no static magnetic screening 
analogous to Debye screening in the electric case, but only a weaker dynamical screening [|J . 
In many cases, this dynamical screening is sufficient to remove the IR divergences |J, || but 
it has been known for some time that it cannot solve easily the IR problem of the damping 
rate [7], at least for non zero T. In a recent paper ||, Blaizot and Iancu were nevertheless 
able to derive a finite result in the T ^ 0, \i = case, by using a non perturbative approach 
to resum the leading divergencies. However they also discovered that the decay law is no 
longer exponential. 

In this Letter, we address the problem of computing the damping rate of quasiparticles 
in degenerate UR plasmas. For the sake of definiteness, we treat the case of a QED plasma, 
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but our results may be trivially extended to the QCD case by substituting to the QED 
coupling e the QCD coupling g, and by taking into account some color group factors. In 
this computation, the basic physical idea is that the collisions of the charged quasiparticle 
with the particles in the plasma are governed by photon exchange, and that one must take 
into account the fact that the photon propagator is dressed by the interactions. Actually, 
this approach is a particular case of the resummation method proposed by Braaten and 
Pisarski ||, which relies on the properties of the so-called "hard thermal loops" [f|, flO 



or "hard dense loops" in the degenerate case [|TTJ. Braaten and Pisarski pointed out the 
importance of a hierarchy of scales, based on the existence of a "hard scale" of order T 
(or /i), and a "soft scale" of order eT (or e/z), with e C 1. When soft scales are involved, 
one must use dressed (or resummed) propagators and vertices instead of the bare ones in a 
perturbative expansion. An important feature of the resummation method is that it leads 
to gauge independent results, due to the gauge independence of the hard thermal (or dense) 
loops. 

Our main result is that in the case T = 0, fi ^ 0, dynamical screening is able to cure 
the IR divergences of the damping rate due to magnetic photon exchange in UR plasmas; 
however, in contrast to the NR case, the damping rate is dominated by magnetic exchange 
and is proportional to (E — /i) , where E is the relativistic energy of the quasiparticle, while 
electric photon exchange gives a contribution proportional to (E — /x) 2 , as in the NR case, 
which may be in fact obtained as a low velocity limit of the relativistic one [12]. Note, 
however, that by convention energies and chemical potentials differ by the rest mass of the 
particle in the NR and relativistic cases. Note also that we use a system of units where 
h = c = ks = 1, and that we follow closely the notations of Jl|. 

Let us now proceed to the derivation of our result. We assume that the quasiparticle 
energy E is hard (this is automatically ensured in the case of a degenerate plasma). The 



damping rate *y(E) is given by the imaginary part of the quasiparticle self-energy S [14 
more precisely 

7 (£) = -J-Tr [Im£(p + i77P)(f + m)] , (1) 

where m is the electron mass, E = (p 2 + m 2 ) 1 / 2 , 77 — > + , and we have used the by now 
standard notation: = (p , p); the lowest order graph for £ is drawn in Fig. ffl. We are 
mainly interested in the contribution of soft photons, so that the electron-photon vertex 
and the electron propagator may be replaced by the bare ones n, [Rf: only the photon 
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propagator need be dressed. 

We perfom the calculation of £ in the imaginary time formalism; then the (free) electron 
propagator is given by 



M^.k)=r^ ( f +m) ^ ) . (2) 

J-oo I'K Kq — %UJ n — /! 

with 

p f {K) = 2ne(k )5(k 2 - E 2 k ) . (3) 

In (0), u n = ir{2n + 1)T is a fermionic Matsubara frequency and e(k ) = k /\k \. The 
(resummed) photon propagator A M ^(Q) is written in the Coulomb gauge 

A^(Q) = 5 M o^o A L (Q) + (<f - - g i g i )A T (Q) , (4) 

where the spectral representations of A^ and Al read 

A /• \ f 00 dq p L {q ,q) 1 

A L z^ s ,g = / ; , 5a 

' 2tt q - iuj s q 2 



oo 
oo 



A T (tu s ,q)= - : — . (5b) 

J-oo Ait go — iuJ s 

In (H) and (|5|), = qj/|q| and u s = 2irsT is a bosonic Matsubara frequency. The explicit 
expressions of the longitudinal and transverse spectral functions pi and px are found by 
taking the imaginary parts of A^ and A^ 

p L (q , q) = 2 Im A L (g + i V , q)= 2 Im — — — ^— — , (6a) 



p T (g , g) = 2 Im A T (g + zt?, g)= 2 Im ^ N , ( 61 > ) 

(go + irj) ~ Q I w p ( x + 2 — m fr~r 



2 ^ v 2 a- 1 , 

where x = (go + irj) / g, ujp = Mj \/3 is the plasma frequency which is related to the Debye 
mass M given by 

The diagram in Fig. p] is now evaluated in the imaginary time formalism (P = (iuj nj p)) 

^ r d^^ 

S(P)=e 2 T^] y — — 7 jU 5' / (z(a; n -a; s ),p-q)7 i/ A^(za; s ,q) . (8) 

The sum over Matsubara frequencies is easily performed when one plugs in ([5]) the 
spectral representations (0) and (||) of the propagators. Taking the imaginary part of E 
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after the analytical continuation iu n + // — > po + if] to Minkowski space, and taking the trace 
in (|l|), one finds for the damping rate, with k = p — q, 

vre 2 f d 3 g r«> dk /•*> dg 

7(E) = 7 (2^ 7-oc ^T P/( " 0) loo (9) 
x (1 + n(gb) - n(fco)) <J(S - fco - go) {[Mo + p • k + m 2 ] 
x Pl(<7o, 9) + 2 [Po&o - (p ■ q)(k ■ q) - m 2 ]p T (g , g)} . 

In (HD, n and n are Bose-Einstein and Fermi-Dirac distribution functions {(3 = 1/T) 

<^ = ^TTT ' ^ = ^(fco-l) + 1 ' (10) 



Eq. @ could have also been derived from kinetic theory |15|], using standard identities 



between the distribution functions fllQl) . The case k < (see @) corresponds in kinetic 
theory to e + — e~ annihilation, which is not IR singular and is even absent in the T = case. 
We thus concentrate on the k > case, which corresponds in kinetic theory to e~ — e~ 
scattering. From now on, we shall also restrict ourselves to the m = and T = case, 



leaving the general case to a forthcoming publication by one of the authors ||12|| . In the 
T = limit, (1 + n(qo)) = Q(qo) and n(k ) = 6(/i — E + q ), where is the step function. 
The qo integration is then limited by 

0<g <£-/u. (11) 

This is also easily seen in kinetic theory, since, due to Pauli blocking, the quasiparticle 
can only scatter into states with energy -E| p - q | such that ^ip-qi < E and furthermore the 
particle on which it scatters must leave the Fermi sea, so that -Ej p - q | > p. Note also that 
the exchanged photon must be space-like: Q 2 < 0, so that the pole part of pl,t does not 
contribute. 

Now, the IR singular contribution comes from small values of the photon momentum q; 
in order to isolate this kinematical region, we follow Braaten ans Yuan |Tj| and introduce 



an intermediate cut-off q* such that e/i <C g* < fi. The "soft" region is defined by q < q*, 
the "hard" one by q > q*: in this latter region we may take the M 2 = limit in the 
denominators of the spectral functions pl,t in (0) M- Let us concentrate on the soft region, 
where we can make the approximation 

£| p _ q , =E-q ~E-p-q. (12) 
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Keeping only the leading terms in (|9]), we find the contribution from the soft region to 
1(E) 



e 



d 3 g 



%oit (E) ^ y / 7^3 ( (<?o) -Q(ji-E + g )) 0(g* ~ q) 

x {p L (qo, ?) + (!- cos 2 9)p T (q , q)} , (13) 

with g = p q = gcos#. It is convenient to use as integration variables go an d g, the 
integration domain D being 

D: {0<q <E-fi; q <q<q*} . (14) 

Then flUf ) becomes (x = qo/q) 

e 2 M 2 f f g 



^soft(E) ~ — / dq dq 

47T JD 



12 . Af 4 7T 2 X 2 



2[q 2 + M 2 Q 1 (x)] z + , 
[2g2 + M2Q 2 (x)] z + ^ ' 



+ To .9 , >fO^ /-_.m2 , i^4 7r 2 a .2 f > (15) 



where 



g 1 ( a? ) = l-^lnl±£, Q 2 (x) = -QW + —L^ . (16) 
2 1 — x 1 — x z 

Note that in the absence of screening (namely, by setting M = in the denominators 
of fll5D ), one would get IR divergent integrals. In general, the integrals in (|TJ) must be 
computed numerically. Fortunately, it is possible to derive an accurate analytical result 
in the physically interesting case (E — fi) <C M. Indeed, it is easy to check that in this 



region one may expand the denominators in (15) in powers of go- Keeping only the leading 



terms, the first denominator in fll5|), corresponding to longitudinal photon exchange, may 
be replaced by (g 2 + M 2 ), which leads to Debye screening. The second denominator in (|T5*D. 
corresponding to transverse photon exchange, may be replaced by 

AqA + fM^ + 8 M 2 q 2 x 2 . (17) 



It can be shown that the last term in ( |TTP gives a subdominant contribution [TJ[], while the 
second term leads to the usual form of dynamical screening J4j - ||. Computing separately 
the longitudinal and transverse contributions, we find, with u* = (q*/M) 2 , 
7 l m . #(E-tf K du ^e 2 M> ( tt 1_\ 

(E -») + —— (E-p) 2 -— - , (19) 



24tt v ' ' 32tt v ' ' V 3g 



where we have only kept the leading terms in (E — //) and 1/q*. 

The total contribution of the soft region to the decay rate is obtained by adding the 
longitudinal and transverse contributions to get 

e 2 e 2 M 2 ( 7T 1 \ 

fso/n 7 24tt V P; 32tt V P; \2M 3 q* 3 ) K ' 

The transverse contribution dominates over the longitudinal one for small values of (E — //) . 

We finally evaluate the contribution from the hard region. Since we are only interested 
in extracting the leading dependence in the fermionic energy of the decay rate we will 
use a simple approach to compute the hard contribution. It is possible to recover bare or 
unresummed perturbation theory to order e 4 by using the spectral densities (|6]) neglecting 
M 2 in the denominators [|]]. This is only valid in the momentum transfer region q > q*. 
Therefore one finds for the hard contribution to the decay rate 



lhard(E) 




After an straightforward computation one finds 

7 ^>^ ( ^(-L-£), (22) 

where q ma x — H is the maximum momentum transfer that it is allowed by kinematics. 

The total decay rate is found just by adding the soft and hard contributions. Then one 
finds that the dependence on the scale q* cancels, as it should. The result is 

1 (E)~—{E-ri + e —{E-tf{—- -L ) . (23) 

In conclusion, we have been able to compute the damping rate of a quasiparticle in 
a degenerate ultrarelativistic plasma, when the fermion energy E is just above the Fermi 
energy /z. This damping rate is dominated by transverse photon (or gluon) exchange and 
proportional to (E—fi). This behavior arises from the combined effect of dynamical screening 
and phase space restrictions due to Pauli blocking. The lifetime r is related to 7 by r ~ I/7, 
and therefore the lifetime becomes infinite as the fermion energy approaches the Fermi 
energy, so that the Fermi sea is stable. 
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FIGURES 



FIG. 1. Resummed one- loop self-energy of the fermion. 



9 



